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We presents a fully quantal version of the Langevin model for the total rate of exoergic ion- 
molecule reactions or inelastic processes. The model, which is derived from a rigorous multichannel 
quantum-defect formulation of bimolecular processes, agrees with the classical Langevin model at 
sufficiently high temperatures. It also gives the first analytic description of ion-molecule reactions 
and inelastic processes in the ultracold regime where the quantum nature of the relative motion 
between the reactants becomes important. 
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The classical Langevin model for exoergic ion-molecule 
reactions 0, Q is one of the most fundamental and pow- 
erful results in the theory of reactions. It has been shown 
to be applicable to a variety of systems and over a wide 
range of temperatures 0, Q • As a model based on the 
long-range interaction, it can be expected to be more ac- 
curate the lower the temperature |4j, until one reaches 
a regime where quantum effects, more specifically the 
quantum effects associated with the relative motion of 
the reactants, become important. 

An experimental study of reactions in this tempera- 
ture regime, often referred as the ultracold regime, has 
recently been realized for neutral-neutral reactions in a 
landmark experiment by the JILA group [Bj]. While it is 
not yet realized for charge-neutral systems, the growing 
ability of making cold molecular samples [6| and manip- 
ulating cold ions implies that it may soon become a 
reality. A sample process, likely among the first investi- 
gated, would include the type 



A+ + B 2 -> AB + B A 
-> A + B+ , 
-> (AB)+ + 



B 



(la) 
(lb) 
(lc) 



where all reactions can be expected to be exoergic if the 
ionization potential of atom A is considerably greater 
than that of atom B. They can proceed with substan- 
tial rates even in the limit of zero temperature as ion- 
molecule interactions are generally expected to be barri- 
erless at the short range [8(. 

The theory presented here, to be called the quantum 
Langevin (QL) model for the total rate of exoergic ion- 
molecule reactions or inelastic processes, gives the first 
theoretical prediction on where and how the quantum 
effects come into play and how the resulting behavior 
deviates from the classical Langevin model. It is another 
application of a new quantum framework for reactions [9j 
that differs considerably from existing formulations [ldj . 
and is used here to further illustrate the concepts behind 
the theory. 

In a conventional quantum theory of reactions [To| , lit- 
tle can be known about a reaction without a detailed 



knowledge of the potential energy surface (PES), the ac- 
curacy of which is often insufficient for quantitative pred- 
ications, especially at low temperatures. This has been 
true even for ultracold atom-atom and ion-atom interac- 
tions, where, to the best of our knowledge, no ab initio 
PES for alkali metal systems has ever been sufficiently 
accurate to predict the scattering length. All potentials 
had to be modified by incorporating a substantial amount 
of spectroscopic data (see, e.g., Ref. [11]). The same is- 
sue becomes much more severe for multidimensional PES 
in reactions, and will remain so for many years to come. 
This difficulty, coupled with the exponential growth of 
the Hilbert space beyond two-body [12], has limited the 
conventional approach to a few simple systems such as 
D+H2, with little hope for more complex systems. 

The multichannel quantum-defect theory for reactions 
(MQDTR), as outline in Ref. [9j], with important moti- 
vations and ingredients that come before it 
offers a different perspective on reactions and inelastic 
processes. It comes from a much wider assertion that 
much can be known about a quantum system, specifi- 
cally its behavior around a fragmentation threshold, sim- 
ply from the types of long-range interactions among its 
constituents. Whatever not yet known can be character- 
ized by a few energy-insensitive parameters, which can 
be further determined from a few experimental measure- 
ments without any knowledge of the short-range interac- 
tion, or even the strength of the long-range interaction. 
This physical picture, which goes back to the ori gina l 
quantum-defect theory for the Coulomb interaction [17], 
has been well established in recent yea rs for atom-atom 
[lsI . [l9j and ion- atom interactions [20||, and to a lesser 



degree for few-atom [21| and many-atom systems [22] . 
Its realization for reactions |9|, as is further illustrated 
in this work, frees the theory from being held hostage 
by the details of PES, while ready to take advantage of 
them when they are available. It also resolves a con- 
ceptual disparity in existing theories of reactions. While 
many classical models 0, |j] are based on the recognition 
of the importance of the long-range interaction, the same 
physical concept gets lost in the conventional quantum 
formulations [lOj. It is primarily due to this omission 
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that they have missed the universality and the simplicity 
in reactions that have been uncovered in the landmark 
JILA experiment 

The QL model for ion-molecule reactions or inelas- 
tic processes, to be presented here, is a special case 
of a more general QL model @ that is formally appli- 
cable to bimolecular processes with arbitrary —C n /R n 
(n > 2) type of long-range interaction in the entrance 
channel. Other than the exponent n, the most impor- 
tant characteristic of such a potential is its length scale 
fi n = (2fiC n /h 2 ) 1 ^ n ~ 2 \ where /i is the reduced mass in 
the entrance channel. It determines the scale parame- 
ters for other relevant physical observables, such as the 
energy, with a scale of se = (?i 2 /2^)(l//3„) 2 , the temper- 
ature, with a scale of s^/fcgj and the rate of reactions, 
with a scale of sk — ^f^fin/fJ- [S|- For ion-molecule in- 
teractions, n = 4 and C4 = aq 2 /2, corresponding to 
the polarization potential, with a being the average po- 
larizability of the molecule and q the charge of the ion. 
The realization of the QL model for this class of sys- 
tems is made possible by a recent reformulation of the 
quantum-defect theory (QDT) for — 1/i? 4 potential [20j |. 
which gives, in particular, analytic results for the quan- 
tum transmission probability to be used in this work. 

In Ref. Q , we have shown that under the Langevin as- 
sumption, corresponding to the assumption of no reflec- 
tion by the inner potential, the total rate of reactions and 
inelastic processes follows a universal behavior uniquely 
determined by the exponent n characterizing the type of 
long-range interaction in the entrance channel. Different 
systems with the same type of potentials differ from each 
other only in scaling. Specifically, the rate constant for 
the total rate of reactions and inelastic processes can be 
written as 

K(T)=s K K^(T a ), (2) 

where sk is the rate scale defined earlier, and K^ n \T B ) 
is a universal function of the scaled temperature, T s = 
T/(s E /k B ), given by 



K, {n) {T s ) = 



dxx 1/2 e~ x W (n) (T s x) . (3) 



Here W^ n '(e a ) is a scaled total rate before thermal av- 
eraging. It depends on energy only through the scaled 
energy e s = c/se, and has contributions from all partial 
waves 



w(")( £s )=^wr ) ( £s ), 



1=0 



(n) 

where Wy is a scaled partial rate given by 
W ( t n \e s ) = (2l + l)T l c{n \e s )/e\l 2 , 



(4) 



(5) 



in which T^ n \ts) — |*/°*^( e «)| 2 is the quantum transmis- 
sion probability through the long-range potential at the 
scaled energy e s and for partial wave 



Aoi) 




FIG. 1. (color online) The scaled partial rates W/ n '(e 3 ) for 
n = 4, corresponding to —1/R 4 type of interaction in the 
entrance channel. 




FIG. 2. (color online) The scaled total rate W (n) (e s ) for n = 
4, corresponding to —l/R 4 type of interaction in the entrance 
channel. 



For n = 4, corresponding to the ion-molecule inter- 
action, the quantum transmission probability, IT (e s ), 
which is the only quantity required to determine the 
universal rate function in the QL model, can be found 
analytically as a part of the QDT for —l/R 4 potential 
H,|2(J. The result is 



r; 



-c(4) 



2M 2 sl [l - cos(2tiV)] 
1 - 2M 2 j cos(2tw) + M e 4 



(6) 



where v is the characteristic exponent and M €a i is one 
of the universal QDT functions for the —1/i? 4 potential, 
both of which are as given in Ref. [2(J . 

Figure [T] illustrates the resulting scaled partial rates 
W\ n) for n = 4. Fi gure shows the corresponding to- 
tal rate W^ n \ We note that the oscillatory structure in 
the total rate is neither a resonance nor an interference 
phenomenon. It is instead a result of the quantization of 
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FIG. 3. (color online) The universal rate function K.' n ^(T a ) 
for — 1/R 4 type of interaction in the entrance channel (solid 
line). The dash-dot line represents the prediction of the clas- 
sical Langevin model, as given by Eq. (|12[) . The dashed line 
represents the QDT expansion as given by Eq. (fTTj) . 



angular momentum, with contributions from a discrete 
set of partial waves peaking at different energies, as illus- 
trated in Fig. [TJ The V\A n ) is related to the total reactive 
and inelastic cross section, cr ur , by <7 ur = (7r^)s£ 3 (e s ), 
in which Eur^(e s ) = /el^ 2 is the scaled total inelas- 
tic and reactive cross section. We note that the oscilla- 
tory structure in is much less visible in the cross 
section, making it difficult to detect experimentally. 

Figure [3] illustrates the universal rate function K.^ (T B ) 
for n = 4, which is the thermal average of the total rate 
shown in Fig. [5] In the ultracold regime of T s <C 1, sim- 
pler analytic formulas for rates and cross sections can be 
derived using the QDT expansion [231 ] of the transmis- 
sion probability. For small scaled energies e s , it can be 
shown from the QDT for n = 4 jH that 



<T-c(4) 
'1=0 



4« 



(4) 1/2 

; ;=o ts 



(1 



i (4) 

L sl=0 



0(e 



5/2\ 



T-c(4) _ 
'1=1 — 



4a (4) e 3/2 



7^ + °( e I /2 ) 



and 



where 



-z; c > (4) = 4ai 4 )4 +1/2 + o(4 +5/2 ) 



(2Z + l) 2 
[(2Z + 1)!!] 



(7) 



(8) 



(9) 



(10) 



is called the scaled mean scattering lengths for a — 1 / R A 
pot ential, after similar quantities for — l/R e potential 
[23[. We note that such analytic expansions could not 



TABLE I. Sample scale parameters for ion-molecule reac- 
tions and inelastic processes. The Pa = (2/iC , 4/7i 2 ) 1 ' /2 is 
the length scale associated with the polarization potential 
— C4/R 4 , where Ca = aq 2 /2 with a being the average po- 
larizability of the molecule, s^/fes = (h 2 /2^)(1/ Pa) 2 /ks is 
the corresponding temperature scale, sk = ttUPa/ ' ix is the 

rate scale corresponding to Pa- It is given here in units of 

" -1 



System 




a (a.u.) 


13a (a.u.) 


s E /k B (K) 


sk 


D+- i H 2 




5.41^ 


99.7 


8.65 x 10- 3 


1.05 


1 H+- 7 Li 2 




216^ 


608 


2.49 x 10~ 4 


6.83 


7 Li+- 87 Rb 2 




553^ 


2610 


1.90 x 10" 6 


4.08 


138 Ba+- 87 Rb 2 


553^ 


8800 


1.45 x 10" 8 


1.21 


40 Ca+- 133 Cs 2 


675^ 


6540 


5.83 x 10" 8 


1.99 


40 Ca+- 40 K 87 Rb 


526^ 


5400 


9.78 x 10" 8 


1.88 


a From Ref. 
b From Ref. 


27]. 
28]. 



have been derived from either the WKB theory [2J], 
or the top-of-barrier analysis [25|, as to be discussed in 
more details elsewhere (26|. Substituting Eqs. (H])-© into 
Eqs. ©-([5]) gives the expansion of K^(T S ) in the ultra- 
cold regime of T s -C 1 , 



JC^(T S ) = Aa^ =Q - lb ^ a po) T i/2 



+18 



(«ifio) 3 + «Hii T s + 0(T*). (11) 



;( 4 ) 



Here a^L = 1 and = 1/225 are the scaled mean 

scattering lengths for I — and 1 = 1, respectively. A 
comparison of this QDT expansion with the exact result 
is shown in Fig. [3] 

At high temperatures as characterized by T s 1, it is 
straightforward to show, from the semiclassical limit of 
the transmission probabilities [lfjj ]. that 



= (4) 



K (4) CT.) 



(12) 



in agreement with the classical Langevin model P, 0]. 
Figure [3] shows that the transition from quantum to semi- 
classical behavior occurs over a temperature range of 
ss/kB, which we generally refers as the van der Waals 
temperature scale. The rate goes from Ask at the thresh- 
old to approximately the Langevin rate of 2sk beyond 
SE/kB- Figure[3]also shows that the oscillatory structure 
present in W^ n ' has mostly been washed out by thermal 
averaging. 

All scaled results can be put on absolute scales us- 
ing a single parameter, the average polarizability of the 
molecule a. It determines both the rate scale sk and the 
temperature scale SE/ks- Sample scale parameters are 
given in Table Q] They are chosen to illustrate that the 
meaning of the ultracold regime, if defined as the range 
of temperatures over which the quantum effects are im- 
portant, is very different for different systems. It covers a 
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much broader temperature range for lighter systems than 
for heavier ones. They are also chosen to imply that we 
expect the QL model to be applicable not only to non- 
polar molecules, but also to small polar molecules such 
as KRb 5] , with the main difference bein g th at its upper 
range of applicability will be more limited [29( . It is inter- 
esting to note from the table that despite wide variations 
of temperature scales for different systems, the rate scales 
are of the same order of magnitude, ~ 10~ 9 cm 3 s _1 , 
which is roughly 100 times greater than those for neutral- 
neutral reactions [9J . 

In conclusion, we have presented a fully quantal ver- 
sion of the Langevin model for exoergic charge-neutral 
reactions and inelastic processes. It is a universal model 
in which different systems differ only in scaling, further 
illustrating the concept that even in a purely quantum 
theory, there are important aspects of reactions that can 
be understood without detailed knowledge of the PES at 
the short range. Such aspects include not only the to- 
tal rate of reactions and inelastic processes, presented 
here, but also the elastic cross section and the total 
cross section, to be presented elsewhere. For a state- 
to-state partial cross section, it can be shown from the 
underlying MQDTR Q that while its absolute value re- 
quires the short-range PES, its energy dependence can 
still be parametrized using the same universal transmis- 
sion probabilities T^ n \e s )- We point out that the QL 
model is applicable not only to molecules in the ground 
or low-lying states, but also to molecules in vibrationally 
highly excited states and to atoms in selective Rydberg 
states (ones with significant quantum defect). In such ap- 
plications, the theory connects with quantum few-body 
physics [30( , and provides a description of their behavior 
outside of the so-called universal regime, a region that 
has been difficult to treat using other methods because 
of the large number of open channels. In a mathematical 
abstraction with even broader implications, the QL mod- 
els presented here and earlier in Ref. Q represents one 
type of universal behaviors that can emerge whenever the 
number of open channels in a set of coupled channel (or 
close-coupling) equations becomes large. It is our belief 
that uncovering and taking advantage of such universal 
behaviors will be the key to a more systemic understand- 
ing of quantum systems beyond two-body. 

This work was supported by the NSF under the Grant 
No. PHY-0758042. 
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